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ABSTRACT 



Plasma balls and rings emerge as fluid holographic duals of black holes and black rings in the 
hydrodynamic/gravity correspondence for the Scherk-Schwarz AdS system. Recently, plasma balls 
spinning above a critical rotation were found to be unstable against m-lobed perturbations. In the 
phase diagram of stationary solutions the threshold of the instability signals a bifurcation to a new 
phase of non-axisymmetric configurations. We find explicitly this family of solutions and represent 
them in the phase diagram. We discuss the implications of our results for the gravitational system. 
Rotating non-axisymmetric black holes necessarily radiate gravitational waves. We thus emphasize 
that it would be important, albeit possibly out of present reach, to have a better understanding of 
the hydrodynamic description of gravitational waves and of the gravitational interaction between 
two bodies. We also argue that it might well be that a non-axisymmetric m-lobed instability is also 
present in Myers-Perry black holes for rotations below the recently found ultraspinning instability. 
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1 Introduction 

The idea of an hypothetical connection between black holes and fluid balls has been the subject of 
several studies along the decades. Nevertheless, for a long time no formal map between these two 
distinct systems was available. (Although the membrane paradigm, whereby a black hole horizon 
is mimicked by a stretched fluid membrane, provided an interesting analogue model with useful 
applications in astrophysical systems). However, the state of the art changed radically in recent 
years. Motivated by the AdS/CFT duality and by Landau's observation that any quantum field 
theory should have an effective hydrodynamic description at high energy densities (where the mean 
free path of the theory is small) , a rigorous duality between gravity and hydrodynamics was finally 
found. Emerging with the path-breaking studies [Il[2], themselves motivated by the ideas of [3]-[H], 
this duality has been since then successfully explored in a series of works [9] -[29]; for a nice review 
and references see [30] . 

The main idea is as follows [H [21 [30] . Start with a black hole solution of Einstein- AdS gravity and 
divide its horizon in several patches. Consider then the approximation where in the tubewise region 
that goes from each horizon patch up to the AdS boundary the geometry is approximately described 
by the line element of a boosted planar black brane. This is parametrized by the temperature T 
and boost parameters /3j along the timelike and boundary spacelike coordinates x*. When moving 
along the several patches one allows the temperature and boosts to vary slowly with the boundary 
coordinates. Obviously such a patched geometry will generically not be a solution of Einstein- AdS 
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gravity. However, it will be a solution at any order of a perturbative expansion of the Einstein- 
AdS field equations as long as at each order n the holographic stress tensor T^^"^ satisfies certain 
equations that turn out to be simply Vi^T/"^ = 0, i.e., it must be conserved (the expansion is done 
assuming that the lengthscales of the temperature and boost variations are much larger than the 
thermal lengthscale of the system) . This stress tensor is found by computing the extrinsic curvature 
of a constant radial surface in the limit where this surface approaches the AdS boundary. The idea 
that the stress tensor must be conserved is certainly not new and was realized in previous studies 
of holographic descriptions of gravity as a quantum field theory on the boundary (namely in the 
AdS/CFT duality). However, with the exception of a few explicit examples, this stress tensor was 
a priori an arbitrary symmetric tensor in d dimensions with d{d + l)/2 components (eventually 
d{d+ l)/2 — 1 if the quantum field theory on the boundary is conformal and thus traceless). These 
are too many degrees of freedom for the conservation equations to be closed. What is new with 
the boosted black brane patching and the perturbative analysis that follows is that in the process 
the form of the stress tensor is highly and uniquely constrained. For example, at leading order 
in the expansion it must have a perfect fiuid form. This reduces strongly the degrees of freedom 
and the conservation equations now constitute a closed system: we have as many equations as 
unknowns. Moreover, there is also a one-to-one map between the bulk gravity solution and the 
boundary system. Working out the conservation equations we further find that they boil down to 
the equations of fluid dynamics. 

These ideas were first developed to get a hydrodynamic description of AdS/CFT duality. How- 
ever they naturally extend to other gravity/gauge theory dualities. Here we will be interested in 
dual systems that have a confinement / deconfinement phase transition [311 [32] . A particular exam- 
ple, that is the simplest one with these properties, is the Scherk-Schwarz (SS) gravity/gauge system 
[531 [M]. The reason to study such a system is (at least) two-folded. First, from a quantum field 
theory perspective, the SS gauge theory is non-supersymmetric, non-conformal and has the above 
mentioned phase transition. These are important properties that it shares with the real world QCD. 
Therefore the SS system can be a useful toy model to infer properties of QCD. In particular, having 
the gravity dual of the SS gauge theory we can use the fact that this is a weak/strong coupling du- 
ality to use the weak field gravity results to learn about the strong field (non-perturbative) regime 
of the field theory. Second, from a gravitational perspective, previous studies of the hydrodynamic 
SS system revealed an interesting but a priori unexpected feature [3l]-[39][l9]. SS AdS gravity 
differs from the globally AdS gravity (some of the main differences will be discussed in Section 
[3]). However, results from SS hydrodynamics are telling us that black holes in this background 
have surprising similarities with globally AdS black holes or even asymptotically fiat black holes. 
We can therefore use the SS hydrodynamic results to learn about these latter black holes in more 
conventional backgrounds. 

Concretely, in the long wavelength regime, the SS compactification of a {d + l)-dimensional 
CFT has an effective description as a d-dimensional fluid dynamics with an equation of state that 
describes the SS plasma. This system is dual to the SS compactification of AdSrf+2 gravity. In 
short, the SS theory has both a confined and a deconfined phase. In the gravity side of the duality 
the deconfined phase corresponds to a black hole solution localized in the infrared region of the 
holographic direction, while the confined phase maps to the AdS soliton solution. They compete 
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with each other for the minimal free energy. In the neighborhood of the confinement temperature 
the two phases co-exist in different regions of the spacetime and are separated by a domain wah 
whose tension provides the Israel junction conditions between the two solutions. In the holographic 
boundary where the fluid lives, the deconfined phase is described by a plasma lump immersed in the 
confined sea phase. The plasma lump has a surface tension at its boundary that is in correspondence 
with the domain wall tension in the bulk. For a more detailed description of this system we refer 
the reader to (32tl38]. 

For the present discussion, we should keep in mind that plasma lumps in d-dimensions cor- 
respond to SS AdSd+2 black objects. For example, axisymmetric plasma balls and plasma rings 
correspond, respectively, to rotating black holes and black rings in SS AdSrf+2 [22] • In a previous 
study, rotating plasma balls were found to be unstable against m-lobed perturbations [38]. Starting 
with a static plasma ball and increasing its rotation, the axisymmetric rotating plasma ball becomes 
unstable first to a 2-lobed ("peanut-like") perturbation and then to m-lobed perturbations with 
m > 2. This conclusion followed from analyzing perturbations of the hydrodynamic equations. 
Usually, such an instability signals a bifurcation point to a new branch of stationary solutions (that 
obeys the symmetries of the instability) in the phase diagram of stationary solutions. Here, we will 
find that this is indeed the case and we will construct explicitly the branch of 2-lobed plasma balls 
in the phase diagram of solutions. 

The paper is organized as follows. Section [2] reviews the relativistic hydrodynamics of (2 -|- 1)- 
dimensional, stationary plasma configurations and displays the foundational equations for the re- 
mainder of the article. The main new results are contained in Section[3l where the non-axisymmetric 
profile of the m-lobed plasmas is found and the phase diagram for stationary plasmas is determined. 
For completeness, we briefly review the previuosly known results for axisymmetric solutions, both 
plasma balls and rings. In addition, we present a simple derivation that precludes the existence of 
m-lobed plasma rings as perturbations of axisymmetric plasma rings. We conclude in Section H] with 
a discussion and some speculations on the gravitational duals of these plasma lumps. Appendix El 
contains an analysis of the non-relativistic limit. 

2 Hydrodynamics of stationary plasma lumps 

In this section we review the relativistic hydrodynamic equations that govern a Scherk-Schwarz 
plasma in a 3d Minkowski background. We will be interested in plasma configurations in mechanical 
and thermodynamic equilibrium. This includes non-axisymmetric configurations. We discuss the 
conserved charges for stationary plasmas that allow to represent the branches of stationary plasmas 
in a phase diagram. We follow closely [36^ I32j. 

2.1 Relativistic hydrodynamic equations 

SS hydrodynamics is an effective description at long distances of SS gauge theory, valid when the 
fluid variables vary on distance scales that are large compared to the mean free path of the theory, 
'mfp ~ 1/Tc- Expanding the holographic stress tensor r^,y in powers of derivatives of u^, one finds 
that at leading order T^j^ has a perfect fluid plus a surface tension boundary contributions, while the 
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dissipation effects due to shear and bulk viscosity plus heat diffusion appear in the next-to-leading 
order term. 

The hydrodynamic equations follow from conservation of the stress tensor, V^T^*^ = 0. This 
yields the relativistic continuity, Navier-Stokes and Young-Laplace equations, 

n^V^p +{p + P)i9 = Ci9' - q^'a^ - V^q^" + 27]a^"'V^u, , (2.1) 
(p + P)a'' = -P^^V^P + C (P^'V^T? + i^^x'^V^n") + 27? (V^a'^'^ - n^a^^V^n,) 

- {qf'Vf.u'' + ^q"" + u^V^g"^ - qf'af.u") , with P'"' = g'"' + u^u'^ , (2.2) 
1 M< 

=aK , with K = /i/V^^n^ . (2.3) 

> 



P - C-d + 2ri { ^'d + u^'n'^Van^ 



Here, u'^ is the fluid velocity, p, P, o", C, and k are, respectively, the density, pressure, surface 
tension, bulk viscosity, shear viscosity, and thermal conductivity of the fluid. The quantities, a^, 
■d, a^'^ and q^ are respectively the acceleration, expansion, shear tensor, and heat flux. P^'^ is the 
projector onto the hypersurface orthogonal to u^^. The fluid boundary is defined by f{x'^) = 0, it 
has unit spacelike normal = d^f /\df\, and h^^'^ = g^^^ — n^^n^ is the projector onto the boundary. 
K is the trace of its extrinsic curvature, and [<5]> = (3< — Q> is the jump on a quantity Q when 
we cross the boundary from the interior into the exterior of the plasma. In the derivation of (j2.3p . 
the constraint that the fluid velocity must be orthogonal to the boundary normal is used (this 
guarantees that the fluid is confined inside the boundary), 

u^n^ = 0. (2.4) 

A fluid with local entropy density s and local temperature T satisfies the Euler relation and 
the Gibbs-Duhem relation given, respectively, by 

p + P = Ts, dP = sdr, (2.5) 

In the SS dual system, we are interested in the long wavelength limit of a Scherk-Schwarz 
compactification of a 4-dimensional CFT. The 3d (non-conformal) plasma that results from the 
dimensional reduction of the 4d conformal plasma has equation of state [32] 

3 1 

P = ^— ^ P + P=3(P-P0), . = 4a.^^J , r=^^J . (2.6) 

with pq and a constants. This equation of state is valid in or out of equilibrium and is normalized 
such that the vacuum pressure vanishesll] 



2.2 Stationary plasmas 

In this subsection we briefly review some general results derived in [36j which are valid for equilib- 
rium plasmas. For later sections, it is important to emphasize that these results are valid not only 
for axisymmetric but also for non-axisymmetric plasma configurations. 

^The simplest possible plasma configuration - a domain wall separating the fluid from the vacuum - is in equilibrium 
when the plasma pressure vanishes. This occurs for a critical temperature Tc = (po/a)^''* |32| . or equivalently for a 
critical fluid density pc = 4po for which P = 0. 
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Stationary plasma configurations are both in hydrodynamical, u^V^P = 0, and thermodynam- 
ical equifibrium u^V^T = 0. Equivalently, they have ■!? = 0, a^'^ = and = 0. Consequently, 
the velocity field of a stationary plasma must be a linear combination of the background Killing 
vector fields and x the stationarity and rotational Killing vectors), 

n = ^(e + Ox). (2.7) 

In particular this implies that a stationary plasma must be at constant plasma temperature T 
related to the local temperature T by the Lorentz factor 



T 

which is the redshift factor relating measurements done in the laboratory and comoving frames. 
Combining the Euler relation (j2.5p and the Young-Laplace equation (equation (j2.3|) without the 
dissipative terms), we can relate the plasma temperature T to a combination of several magnitudes 
at the fiuid surface, 

„ aK +0 , , 

T= —. (2.9 

In the plasma lump/black hole duality, T maps to the Hawking temperature of the horizon. A 
stationary plasma must also be rigidly rotating, i.e., with constant angular velocity 0,, which in 
turn maps to the horizon angular velocity in the duality. 

Summarizing, the equations that a stationary fiuid must satisfy are given by (j2.1J) - ()2.3p with 
= 0, a^'^ = and = 0. In particular the Young-Laplace equation ()2.3p reduces tqj 

P<-P>=aK. (2.10) 

For the SS plasma in equilibrium, it further follows from ()2.8p and ()2.6p that the pressure and 
energy density satisfy the relations 

P = ^7^-po, P = P*l^ + Po, (2.11) 

where = 3aT^ is a constant for any given plasma configuration. 



2.3 Conserved charges 

The constituent fiuid of the plasma object has local energy density p, pressure P, velocity u^, local 
entropy density s, and local temperature T. These local quantities provide the information we need 
to compute the thermodynamic quantities (energy, angular momentum, entropy, temperature) of 
the stationary configurations (plasma balls, plasma rings and plasma peanuts). 

Here, we restrict our discussion to fiuids that live in a (2+l)-dimensional Minkowski back- 
grouncjfl with stationarity timelike Killing vector ^ = dt and spacelike Killing vector x = {^or 

^The last term on the LHS of (|2.10p vanishes in equihbrium because use of n^n'^a-fj,^, — imphes that u'^n°'V — 
■^We use polar coordinates {t,r,(j}), and the non-vanishing affine connections are F^^ = — r and F^^ = F^^ = 1/r. 
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general results see [36]). We can then foliate the spacetime into constant t hypersufaces T,t and 
is their unit normal vector. Then, given any Killing vector ^p^, one can define the associated 
conserved charges Q[ip] = Jj^dV T^^^'^fp'^, where dV is the induced volume measure on S^. The 
fluid velocity is given by ^7f\ . i.e., u^" = (5''* + QS^'f') with 7 = (l - r"^ n'^y^^'^ . The energy 
and angular momentum of the plasma associated, respectively, with the Killing vectors ^ and x 
then 

E[V]= f [{p + Pm-uf + iC-C)P]dV-a[ eeh^u\dmf)dV, 
Jv JSt 

J[V] = f {p + P){^.u){x-u)dV-a f h^^exldmDdV . (2.12) 

Note that for axisymmetric solutions (plasma balls and plasma rings) the boundary term in J 
proportional to a vanishes. It is however present for non- axisymmetric solutions where x ■ ^ 7^ 
{i.e., when the fluid boundary is not invariant under the action of x) |36j . The total entropy of the 
fluid is the conserved charge associated to the entropy density current su^, 

S = - [ dV{^-u)s= [ dVjs. (2.13) 
Jv Jv 

3 Stationary solutions: balls, rings and lobed plasmas 

In this section we are interested in rigidly rotating equilibrium configurations in a (2+l)-dimensional 
Minkowski background. In addition to the two families of axisymmetric solutions composed of 
plasma balls and plasma rings there exist also non- axisymmetric configurations or m-lobed plasmas, 
where m runs over the set of positive integers excluding 1. The plasma balls and rings were discussed 
in detail in [32] . Because we will later study the phase diagram of these solutions we review here 
their properties, following closely |32j . 



3.1 Axisymmetric solutions: plasma balls and plasma rings 

Consider plasma configurations in a 3d Minkowski background parametrized by coordinates {t, r, (p). 
The axisymmetry requirement demands that the boundaries of the plasma depend only on r. Each 
boundary is thus defined by the condition (a specifies a particular boundary in the case where more 
than one is present) 

/(r)=r-i?„ = 0, (3.1) 

and has unit normal = = 6^r- Its extrinsic curvature is -ftT = Following ^32j, it is 
convenient to frame our discussion in terms of the dimensionless variables. 



Q = ^^ f=P^^ y = Qr = Qf, (3.2) 
Po o- 



and also to use dimensionless thermodynamic quantities. 



£ = S = 4^, f = r(°)'". (3.3) 
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In 3d the possible axisymmetric stationary configurations are the plasma balls and the plasma 
rings: 




• Plasma balls are characterized by having a single axisymmetric outer surface at r = Rq and by 
P> = 0. Its properties can be found using the equation of state (|2.1ip . and the Young-Laplace 
equation (|2.10p . It follows that plasma balls must satisfy the condition (with < Vq < 1) 

{l-vlf. (3.4) 



Plasma rings have an axisymmetric inner surface at r = (where P< = 0) , in addition to the 
outer surface at r = Ro (where P> = 0). The equation of state (|2.1ip . and the Young-Laplace 
equation (|2.10p introduce a constraint equation for these objects, 

(3.5) 

It constrains the three variables Vq, vi and Q as, e.g., Vi = Vi{vo,Q). An inspection of it 
concludes that there is a minimum Vq, call it v*, above which (j3.5p is valid ^2]. So, plasma 
rings exist only for t^o ^ ^o- ^^^t there are two families of black rings. One is called the 
fat plasma ring and exists for Q < Vi < v* (where v* < v* is such that the derivative of 
of the RHS of (13.5P w.r.t. Vi vanishes). The second, dubbed as thin plasma ring, exists for 
v* < Vi < 1. At Vi = V* the two families meet at a regular solution. The condition Vi > Q 
guarantees that the temperature is non-negative. 

The dimensionless energy, angular momentum and entropy of the plasma rings are, respectively, 

<Vo - Vf) - (vt - Vf) + ^^{Vo + ^i) - n{vl + Vf) 




E 
J 



02 

2{vt - vf) + 2n{v^ + vf) 
4 



02 




(3.6) 



and the corresponding charges for the plasma balls can be obtained from these formulas by taking 
Vi = 0. 

3.2 Plasma peanuts and m-lobed configurations 
3.2.1 Boundary equation for lobed plasmas 

According to the discussion in Section 12.21 the stationarity restriction can accommodate non- 
axisymmetric configurations as long as they are in rigid rotation. Taking this into account, we 
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drop the assumption of axisymmetry and define the surface of the plasma lump through the con- 
dition 

f{t,r,cl)) = r-Rit-(P/n) = 0. (3.7) 

The boundary's unit normal is 



n„ = \df\-'{-R'6^ + 6;; + ^6^ 



\0f\ 



(3.8) 



where R' = with x = t — (p/Q. The condition that no fluid flows through the surface, 

namely u'^n^ = 0, is manifestly satisfied by the ansatz (j3.7p when the velocity field takes the form 
u'^ = j{5f +^16^). The trace of the extrinsic curvature of the boundary in the general rigid rotation 
case is 



K = Q- 



-RR" (1 - a'^R^) + (2 - n'^R^) + Q2ij2 

[O^R^ + i?'2 (1 - J]2^2)]3/2 

Introducing the dimensionless quantities (13. 2p and 

ij) = —VLx = 4> — Q.t , Vo{t(;) = Q,R{x) , 'yo(V') = R'{x) , k 



p* 
3po 



the Young-Laplace equation (j2.3p , after employing ()2.6p and (|2.8p , becomes 
vMl-vl)-v',^{2-v^,)-v^, , ^ 

n 



This equation admits a first integral, namely the quantity 



+ ^[kil-vl)-'-l] 



^[k{i-v'^r'-vi]^Q 



(3.9) 
(3.10) 

(3.11) 
(3.12) 



is a constant independent of ip. 

Now consider the stress tensor for a perfect fiuid with surface tension a. The relevant T^jy 
components to compute the charges are 

a n'^R^ + R'^ 



Tu = i' (p + i^Vp)e(-/) + 



t<j> 



Y^r' (p + p)Q(-f)-a 



^R y/ww+wnr^ww) 

RR'^ 



(3.13) 



^1^2^2 + ^/2(1 _f^2^2) 

The induced measure on a constant t hypersurface used to compute the conserved charges in 



(121^ and (l233|) is dV 



§~^vdvdip. Using (j2.12p and (j2.13p . the energy, angular 



momentum and entropy of the lobed plasma are then 



E 



J 



S 



1 



2lT 



27r02 
1 



dxp 



1,2 _|_ y'2 

I/q 



2n 



ttO,'^ Jo 
2 



dip 



(1_^2)^ ■ V^2+^.2(l_^2) 
kvf '^-2„/2 



+ 



Jl_^2)2 ^^2 + „/2(l_„2) 
2- fc3/4„2 



1-t;^ 



(3.14) 
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where Vo = Vo{tp) is the velocity of the plasma peanut boundary that solves the Young-Laplace 
equation and v'^ = Note that by setting Vo{ip) = 0, and k = (l + Q/vo^ (l - (see (f3^ ) 
we get the plasma ball charges (equation ()3.6p with Vi = 0). 

Lobed plasma configuration are those whose boundary profiles satisfy equation (|3.1ip (or (|3.12p ). 
Next we find these solutions in a perturbative expansion around the plasma ball bifurcation point. 

3.2.2 Lobed plasmas from perturbed analysis of plasma balls 

As a warm-up exercise that will also be very useful to check our numerical results against, we find 
lobed plasma configurations by performing a perturbative analysis of the plasma balls around the 
bifurcation point to the new phase. 

Consider a small non-axisymmetric perturbation of the plasma ball, 

Q = n + euj + 0{e'^), k = k + ek + 0{e^) . (3.15) 

where the hat refers to the (constant) unperturbed quantities, e ^ 1 is the expansion parameter, 
and the unperturbed plasma ball radius Vq necessarily obeys p.4p . which we write here as 



k{l-vt)-'-l . (3.16) 



Note that a consistent perturbation requires that we disturb not only Vq but also and k, as done 
in (|3.15p . Linearizing the Young-Laplace equation (|3.1ip we then get 

v" + m'^u = /\, (3.17) 

with 



2 1 4:kvl A - 

= ^ + 13^7 > A = Wo 



k-{l-vl) 



2^2 



17 K 



(3.18) 



I-Vl J](l-l?2)4' 02(i_^2)3 

The solution for the perturbation is oscillatory, 

z/(V') = ^sin(mV' + /3) + A/m^ , (3.19) 

and imposing periodicity ip ^ ijj + 2Tr implies that m ^ 7L. The plasma peanuts correspond to 
m = 2. The definition of the small expansion parameter, 

e = ^ [ # [voW - Vo] z^(V') , (3.20) 



^'o JO 

provides a normalization for the function z^(V'): 

l2 



d^[i/(^)]^ = 1. (3.21) 
This in turn may be used to fix the amplitude A as 







vr 
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However, A does not enter the first order expansions of tlie various tliermodynamic quantities so 
we shall not need it. 

Note that combining eqs. (j3.16p and (|3.18p we obtain the constraint 



m 



(3.23) 



Recall that for a given energy, the unperturbed plasma ball quantities Vo and Q are related by the 
first equation in (13. 6p . Relation (j3.23p then determines the value of Vq (or 0) where the m-lobed 
plasma bifurcates from the plasma ball branch for a given energy. In the non-relativistic limit, 
^ 1 and -^^^ <C 1, eq. ()3.23p reduces to AR^Q'^po/a = rr? — 1 and agrees with the bifurcation 
point found in [38jEI 

Having identified the bifurcation point we can now determine the slope of the new phase branch 
curve in the neighborhood of the bifurcation point. Expanding to first order the dimensionless 
energy, angular momentum and entropy of the lobed plasmas (given by eq. (j3.14p ) around the 
plasma balls, we get 



E 



J 



+0(8^) , 

Jpb + e 



(1 



2 

K — — 








00 + 


[(■4) 







- + 



(1-1}2)3 



kv. 



^3(1 



Spb + e 



r?2(l 



2k Quj 8A 

k h m?{l - 

3k 8a; 8A 

k h m?{l-x 



(3.24) 



Here, the subscript pb is appended to quantities pertaining to plasma balls (see eq. (13. 6p for their 
definitions) . 

In the next subsection we shall obtain the phase diagram at fixed energy. Since we focus 
on 2-lobed perturbations of plasma balls, we require that the energy E of the non-axisymmetric 
configuration actually equals Epi,. Thus, the energy constraint reads 



rr2\2' 



2 



1 + 



+ 



k{3 



(1 



-:2\2 







00 + 


[(■4) 







+ 



2A 



0. 



(3.25) 



Together with (j3.18p . the energy constraint yields a linear relation between the first-order pertur- 
bation quantities. 



00 ■ 



(1 



:v2\2 



(9 



39^2 + 549^ 



27v^ + iv^ 



Vo{-3 + Wv"^ - 8vt + 2v^) 



(3.26) 



''There is an interesting feature associated with this value. In [38] the classical critical rotation was found studying 
perturbations of the dissipative hydrodynamic equations H2.1|l - (|2.3|l in their non-relativistic limit. Curiously, if we 
repeat the perturbation analysis starting this time with the non-dissipative fluid equations we find that the value of 
the critical rotation shifts to 4_Rof2^po/o" = m{m+ 1). So, even a vanishingly small viscosity has a dramatic effect on 
the critical rotation at which instability sets in. Viscosity introduces a singular-limit, where the "limit of the theory 
is not the theory of the limit". In classical fluids this property is very familiar and well studied [4UI I41j . 
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where we have used eqs. ()3.16p and ()3.23p to ehminate k and in favor of (and replaced m = 2). 

At this point we have not yet fully determined the first-order quantities {iv('0), cj, k}, but to 
compute the slope at which the plasmapeanut branches off from the plasma ball curve in the S — J 
phase diagram we need not go furtheio. The reason for this is that the slope (at the bifurcation 
point) is given by 

(~\ peanut / dS \ 

(3.27) 

and therefore the only terms for S and J in eq. (j3.24p that contribute are those linear in e. Taking 
the ratio, the undetermined first-order quantity that was left cancels out, leaving an expression as 
a function of only. The final result for the slope of the plasma peanut at the bifurcation point 
in the S{J) phase diagram is 




peanut „ 

" " (3 - Avlf' [(1 - ^T2)(3 - 11^72 + 4^4)]3/4 • (3-28) 

For E = 40, eqs. ()3.6p and ()3.23p determine the bifurcation point to occur at Vo ~ 0.378 {VL ~ 0.14). 
Plugging this value into eq. (j3.28p we obtain dS/dJ ~ —0.142. This perturbative result will be 
compared with the numerical result of the next subsection and a good agreement will be found. 



3.2.3 Phase diagram of stationary plasmas 

We have integrated numerically equation (j3.1ip to find lobed-like solutions. In the numerical 
search, we fix the energy at ii^ = 40 it 10^'' (for simplicity, since it allows comparison with previous 
investigations |32l 1351 138] ), and vary Q. There are three parameters to be varied in the numerical 
integration: Q, k and Vo{0). For each 0, we search for an array of k,Vo{0), which yield m-lobed 
solutions with ?;o(0)/?;o(27r) = 1 it 10"''. Among those, we select the one for which ii^ = 40 ± 10""^. 
We started the integration close to the plasma ball configuration, by looking for a two-lobed solution 
with i'o(O) ^ Vo{Tr). We find that the bifurcation to plasma peanuts occurs at the critical values 
J ~ 22.2, Q ~ 0.16, S ~ 32.7, and T ~ 1.0. Our results are summarized in Figures [T]l2j 

In Figure [H we show typical two-lobed configurations for different values of rotation parameter, 
Q ~ 0.13,0.12,0.09 respectively. As $7 decreases, the configurations get more and more distorted, 
and for small rotations they assume the form of peanut-like shapes, as in the rightmost panel of 
Figure [TJ Once one knows the configuration shape it is straightforward to compute the energy, 
entropy and angular momentum from eqs. (j3.24p and to establish the phase diagrams S{J), il.{J) 
and T(J). These diagrams are shown in Figure [H where we include also plasma balls and plasma 
rings [32]. For the particular energy we are considering, the 2-lobed configurations have a slope of 
about dS/dJ ~ —0.149. This is in very good agreement (~ 5%) with the perturbative construc- 
tion prediction dS/dJ ~ —0.142 of the last subsection. We have constructed numerically other 
configurations, including classical configurations, for which the velocities are very small, and we 
have obtained good agreement for the slopes predicted in [38] . 

^In the non-relativistic case the perturbative analysis must be extended to second order |38) . 
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Figure 1: Different configurations of a two-lobed plasma for three different values of angular velocity 
Q = 0.13,0.12,0.09. The energy is kept fixed at £" = 40 and the bifurcation rotation occurs at 
~ 0.16. (The rotation axis is normal to the page). 
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Figure 2: Phase diagrams S{J), il(J) and T(J) with the two-lobed, "peanut-like" configurations, 
along with the plasma balls, and thin and fat plasma rings found in [32j. The energy is fixed at 
^ = 40 as in [32] . 

One particularity of the plasma peanuts (common to the thin rings) may be noted here: for 
these configurations an increase of angular momentum produces a decrease of angular velocity, as is 
evident from Figure [2j Such a phenomenon is made possible because the moment of inertia grows 
accordingly and this is apparent in Figure [TJ 

3.3 No m-lobed plasma rings 

A natural question that emerges from the plasma ball analysis done so far is whether or not it is 
possible to have also m-lobed plasma rings. In a previous study this question was addressed by 
perturbing hydrodynamic equations (|2.ip - (|2.3|) with an ansatz that satisfied the symmetries of the 
hypothetical m-lobed plasma rings ^38j. The analysis concluded that plasma rings are stable against 
this particular sector of perturbations. We now revisit this problem from a different perspective. 
The outcome of the present analysis will reinforce the conclusion of the previous study. 

Motivated by the success of the analysis done in subsection 13.2.21 the idea is to ask if we can get 
m-lobed plasma rings as perturbative solutions around a possible bifurcation point in the plasma 
ring branch. It follows from the Young-Laplace equation (j2.10p and equilibrium equation of state 
(12.11j) that the outer boundary of a hypothetical m-lobed plasma ring has to satisfy (13. lip , and its 
inner boundary must obey a similar relation with the trade Vo Vi for fixed {k, fl}. The fact that 
the parameter k = p*/(3po) must be the same for both boundaries in particular implies that such 
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a solution would have to satisfy the relation, 



k = {i-v'^y\i-n 



[V2 + V'^^l-V2)] 



3/2 



In the axisymmetric case, = = v[, this relation reduces to the equilibrium condition (j3.5p . 

Consider now the (natural) possibility that an m-lobed plasma ring could emerge from a slight 
perturbation of an axisymmetric plasma ring. Following a similar strategy as in section 13.2.21 we 
linearize (j3.29p by perturbing the boundaries as 

Vai^j) = Va[l + ei^aW + 0{e'^)] , with a = o,i, (3.30) 

and Q and k as in (j3.15p . Again, the hat refers to the (constant) unperturbed quantities. To leading 
order, we get (13. 5p and in the next-to-leading order we find the two equations, 

u'^ + m\va = , with a = o, i , (3.31) 



and 



2_1 ^^^a A_- 



k-{l-vl)' 



UJ — Qk 



±1. (3.32) 



J- - i2(l-i;^) (1 - t!^) 

Each of these equations admits an oscillatory solution of the form (|3.19p . m-lobed plasma rings 
can exist only if the number of nodes in the outer and inner boundaries is the same, i.e., if the 
condition 

rrio = m-i , mo,miGZ (3.33) 

is satisfied. In this relation and (|3.32p . the unperturbed quantities are constrained as discussed in 
subsection 13.11 e.g., one must have Q < Vi < 1 and k = k{i},Vo) = k{0,,Vi) as indicated in (j3.5p . 
An inspection of p.32p concludes that the equality (|3.33p is verified only when vi > Vq (see Fig. [3|) . 
This clearly violates the basic assumption that the inner radius of the plasma ring must be smaller 
than its outer radius. 

This analysis therefore rules out the possibility that lobed plasma rings could exist as a branch 
of solutions emerging perturbatively from the axisymmetric plasma rings. There is a clear physical 
reason for why the plasma rings are stable against "lobing", in contradiction with plasma balls. 
The additional angular momentum acts to push the fluid radially outward. For the plasma ball 
the only way this can be done while preserving the energy is by breaking axial symmetry, whereas 
plasma rings can accommodate such a change by becoming larger and thinnejfl. 

Reference |38j looked for the possible existence of these solutions following a different path. 
It analyzed the stability of a plasma ring against m-lobed perturbations such that these were in 
phase both in the inner and outer boundary deformations. No unstable mode was found. These 
two negative results following different approaches seem to suggest that the existence of m-lobed 
plasma rings is ruled out. 



^This argument is due to V. Hubcny, to whom we thank. 
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Figure 3: as a function of Va for the three possible regimes of 0. The sohd Hne represents 
m^{^l,Vo), while the dashed line represents mf{i},Vi). The case ruo = mi = 2 is possible only for 
large rotations and when vi > Vq- This violates the plasma ring condition that the inner radius 
must be smaller than the outer radius. 

3.4 The hydrodynamic regime 

Relativistic hydrodynamics provides a good effective description of the deconfined plasma phase of 
J\f = 4: Yang Mills theory compactified down to d = 3 on a Scherk-Schwarz circle only if certain 
conditions are satisfied [32]. First, hydrodynamics is by definition valid when the thermodynamic 
quantities of the fluid vary on a lengthscale large when compared with the mean free path, ^mfp; 
of fluid constituent particles. In our case £mfp ~ ^ good estimate for the valid 

regime is obtained when the maximum fractional rate of change of the fluid local temperature, 
^1 ^ drln^l (recall that T = T^) is much smaller than i~l . This occurs for -pL^o ^ 

I Imax ' Imax ^ mtp 1— I'd 

This condition is satisfied by a wide range of plasma balls and rings as long as we are away from 
extremality [32]. Second, the analysis done so far and onwards assumes that surface tension is 
constant, a = cr{Tc), when in fact it is a function of the fluid temperature at the surface. This 
assumption is valid when T/Tc ~ 1 at the boundary surfaces. This is the case for an extended 
range of energies and angular momenta as long as we do not approach the extremal conflgurations 
|32j too closely. Finally, the boundary of the plasma is treated as a delta-like surface when in fact 
it has a thickness of order T^^. So the analysis is valid when the boundary radius is everywhere 
large when compared with T~^, {Rq, Ri, Rq — Ri} S> This is satisfled if the plasma energy is 
large and again if we are away from the extremal conflgurations [32]. 

4 Gravity duals. Discussion 

Since plasma balls are unstable above a critical rotation rate, the holographic dual SS AdSs black 
holes should also be unstable against m-lobed perturbations. This instability then provides a 
mechanism that dynamically bounds the rotation of these black holes. From the plasma results, we 
expect this bound to be well below the rotation where the extreme solution is reached (for which 
T = 0), and below the minimal rotation where SS AdSs black rings can exist. Very little is known 
about SS AdSs black holes, so we are unable at this point to explicitly check the existence of this 
instability and rotation bound in the gravitational system. 

In the fluid, we explicitly found the new branch of stationary non-axisymmetric plasmas bi- 
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furcating from the plasma ball curve at the point where the instability becomes active. In the 
simplest m = 2 case this is a plasma peanut. In the entropy vs angular momentum diagram (for 
fixed mass) the plasma peanuts have higher entropy than the plasma ball. In fact, for the range of 
angular momenta for which plasma peanut solutions exist, they are entropically dominant over all 
the stationary configurations. 

In the dual gravitational system, m-lobed plasma balls are in correspondence with m-lobed black 
holes. However, the latter can only exist as a long lived object, but never as a stationary solution. 
Indeed, a rotating non-axisymmetric black hole has a quadrupole moment and thus necessarily 
radiates away .ts lobed deformationfl The expectation is that this emission of gravitational waves 
will proceed until the system is back to an axisymmetric configuration spinning at a rate below the 
critical rotation where the instability kicks in. 

This point deserves a few more words. Consider first global AdS(i_|_2 spacetime. Its boundary is 
the Einstein static universe x S'^. This boundary behaves effectively as a reflecting wall or box. 
Therefore, in this background one could eventually have a spinning non-axisymmetric black hole 
surrounded by rotating radiation. The rate at which the black hole would be radiating could in 
principle be balanced by the absorption of previously emitted radiation that is reflected from the 
boundary. An equilibrium solution with a black hole plus orbiting radiation is a priori conceivable 
in this background. On the other hand, SS AdS has boundary X M^-i X Sgg and behaves 
quite differently. In this background, only propagation along the holographic radial direction faces 
reflective boundary conditions. Along the other spacelike directions, namely those parallel to the 
holographic boundary (where the dual fluid lives), we have asymptotically flat boundary conditions. 
It then follows that waves emitted by a black hole along the radial direction will be reflected back 
but radiation emitted along the boundary directions will leak towards infinity. Therefore, the m- 
lobed plasma balls are dual to long-lived lobed SS AdS black holes that decay slowly, i.e., that 
loose angular momentum along the boundary directions and cannot be stationary. 

The existence of the lobed instability and of the stationary plasma peanuts and related families 
exposes the need to explore one point in the fluid/gravity correspondence that is still poorly under- 
stood: what is the fluid description of gravitational interactions and gravitational radiation? At the 
approximation level at which the duality is currently understood, hydrodynamics can only describe 
gravitational systems where these phenomena are suppressed. Indeed, two plasma balls do not 
interact, their collision is not accompanied by radiation emission and a non-axisymmetric plasma 
ball does not radiate. In short, the fluid description only captures the interaction of two black 
holes in the limit when their separation is much larger than the AdS radius, and the gravitational 
emission in the limit where its wavelength is much larger than the AdS scale. 

It is known that the gravitational radiation and interaction corresponds in the plasma to the 
emission of glueballs. It is also known that in the large limit of confining SU {N) gauge theories 
that have a gravitational dual, the energy loss in this emission process is suppressed by a 
factor when compared with the plasma ball energy density. For this reason we know that our plasma 
results describe a necessarily long-lived SS AdS lobed black hole where gravitational emission occurs 
at a very slow rate. It would nevertheless be very interesting to work out more quantitatively the 

static non-axisymmetric black hole was argued to exist in Kaluza-Klein gravity, i.e., that asymptotes to 
M'' X S\ for d > 4 gl]. 
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fluid description of the gravitational emission and interaction processes, to accommodate also in 
the formulation cases where these phenomena are not suppressed. However, it is conceivable that 
such processes might force one to leave the regime of validity of the hydrodynamic description and 
therefore may be hard to captur^. 

Although we have only worked out the lowest dimensional case, higher dimensional SS AdS 
theories also have a fluid description. More specifically, for any d > 3, the SS compactification 
of {d + l)-dimensional CFT has a d-dimensional fluid dynamic description which is dual to a 
SS compactification of AdSrf4.2. An interesting evolution is expected as we climb the dimension 
ladder. As described above, the asymptotics of SS AdS gravity roughly interpolates between AdS 
and Minkowski asymptotics depending on the spatial direction we look at. But as we increase the 
number of dimensions the "AdS-like" radial direction is kept and we are only adding more and 
more "flat-like" boundary directions. Therefore, we antecipate that as d grows the fluid results 
should describe black holes whose properties increasingly resemble those of asymptotically flat 
black holes \36\ I39j . The known results confirm this expectation, as we briefiy review next. Start 
with d = 3. The phase diagram for plasma balls and plasma rings in 3d (see Figure [2]) is similar to 
the phase diagram expected for black holes and black rings in AdSs [32]. In particular, in contrast 
to the asymptotically flat solutions, the black rings have an upper bound on the rotation above 
which they would not fit in the AdS box [32j. So, the 3d plasma results indicate that the SS 
AdSs black objects behave similarly to AdSs black holes. When we increase the dimension by one 
unit two interesting new features emerge [32^ [35] . (The following discussion only highlights some 
properties of interest for our purpose; we ask the reader to see [35] for details). First, the plasma 
rings have now unbounded angular momentum. That is, SS AdSg black rings behave much like 
6d asymptotically flat black rings. So, quite amazingly, in what concerns this particular feature in 
the fluids d = 4 sets already the critical dimension for the above mentioned transition between an 
asymptotically AdS and asymptotically flat behavior. Second, from the gravitational perspective, 
we know that asymptotically flat (singly spinning) black holes in c? > 6 have no bound on their 
rotation. The transition from moderately rotating black holes into ultraspinning ones is set by 
the critical rotation where the so-called ultraspinning instability becomes active |45l |46j!|. In the 
phase diagram of solutions this marks a bifurcation point to a new branch of axisymmetric pinched 
black holes whose horizon is distorted with ripples along the polar direction. Returning to the 
fluid results, in d = 4 a third stationary solution — the pinched plasma balls — is indeed found, in 
addition to the plasma balls and rings [32l [35] . 

We have not worked out the technical details but in d > 3 the m-lobed instability must also be 
present in the plasma balls. The reason being that it is well established that unstable phenomena 

*We thank M. Rangamani for pointing this out to us. 

^In more detail, when rotation keeps increasing in the ultraspinning regime the horizon flattens in the neighborhood 
of the axis poles and the geometry becomes more and more like the geometry of a black membrane [IS]- But the 
latter is unstable against the Gregory-Laflamme instability when their length along the extended directions is bigger 
than their transverse radius [43]. Therefore, at the critical rotation where we enter the ultraspinning regime (and 
well below the limit where the geometry becomes a black membrane) the rotating black holes are expected to become 
unstable against the naturally dubbed ultraspinning instability [35] . A recent study confirmed that this is indeed the 
case [IB]. In the fluid description, the Gregory-Laflamme instability maps into the Rayleigh-Plateau instability that 
is responsible for the pinch-off of long plasma tubes or branes into plasma balls [441 136) . 
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that are already present in classical fluid dynamics cannot cease to exist in relativistic hydrody- 
namics. Now, in [1^ it was found that the lobed instability is indeed active in higher dimensional 
classical fluids. Therefore it should persist in the relativistic regime. The interesting question is 
then whether this instability appears at lower or higher values of rotation than the ultraspinning 
one. We argue that it should appear for lower spins. The argument is simple but seems to be 
robust: in 3d the former instability is present while the latter does not appear. It would become 
active at rotations higher than the upper bound for the angular momentum. When we increase 
d this cap in J disappears and the ultraspinning makes its appearance but, it should do so at a 
rotation larger than the 2-lobed critical rotation. 

As emphasized several times, the plasma results are strictly dual to SS AdS gravity. But as 
described above, in certain regimes, we also have or expect strong similarities between black holes 
in this background and those that are globally AdS or asymptotically flat. It is therefore compelling 
to conjecture the possibility that a non-axisymmetric m-lobed instability might also be present in 
higher dimensional Myers-Perry(— AdS) black holes. We do not find any reason to discard it a 
priori; a numerical investigation of this hypothetic instability might reveal interesting results. 
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Appendix 

A Non-relativistic lobed configurations 

In this Appendix we review the classical hydrodynamic equations that describe non-relativistic 
fiuids and, in particular, lobed configurations |48]-[50]. As a consistency check, we confirm that 
the non-relativistic (NR) limit of our equations for the plasma peanuts reduces to these classical 
equations. 

The classical analysis of rigidly rotating fluids is standardly done in the rotating frame of 
reference. The advantage of this frame is that the (unperturbed) velocity of the fluid vanishes. In 
this frame, the classical Navier-Stokes, continuity and Young-Laplace equations are respectively 
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given by 



atv + (v • V)v = --VP - 2n X V + ([fi X r|) , 
P ^ 

dtp + v ■Vp + pV - v = 0, 

P<-P> = aK, with K = V-n, (A.l) 

where v is the spatial fluid velocity, and V the spatial gradient. In the Navier-Stokes equation, the 
two last terms describe, respectively, the Coriolis and the centrifugal acceleration contributions. 
We are interested in stationary fluids so we did not include dissipation terms that vanish for 
these configurations. We consider incompressible fluids (i.e., with vanishing convective derivative, 
dp/dt = dtp + V • Vp = 0) and thus the continuity equation reduces to the statement that the 
velocity is a solenoid vector, V ■ v = 0. Note that in the classical Young-Laplace equation, K is the 
spatial or mean curvature (and not the spacetime curvature). 

For a fluid in rigid rotation, the Navier-Stokes equation yields for the pressure: P = ^pfl'^r'^ + C , 
where C is a constant. To make later contact with the relativistic result, we choose to replace the 
constant quantities {p,C} by {po,k}, with the relation between them being C = po{k — 1) and 
p = 4:pQk. A stationary fluid with a non-axisymmetric boundary profile described by f{r,(j)) = 
r — R{(f>) = has mean curvature K = '^|j^2^^j;/2)?/^ • Introducing the dimensionless quantities 

i^ = L voW = nRiij), n = — , (A.2) 
n po 

the Young-Laplace equation can then be written as 

'-frr^^ (A.3) 

This equation agrees with the NR limit of ()3.1ip . as it should, and gives the equation for the profile 
of stationary classical lobed configurations. It has the first integral, 

obtained by integration of (IA.3|1 . By curiosity, note that this is an example of a case where "the 
limit of an expression is not the expression of the limit" since (IA.4P is not the NR limit of the first 
integral (l3T2]l . 
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